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ABSTRACT 

We consider the (2n + l)-dimensional euclidean Dirac operator with a 
mass term that looks like a domain wall, recently proposed by Kaplan 
to describe chiral fermions in 2n dimensions. In the continuum case we 
show that the euclidean spectrum contains no bound states with non- 
zero momentum. On the lattice, a bound state spectrum without energy 
gap exists only if m is fine tuned to some special values, and the disper- 
sion relation does not describe a relativistic fermion. In spite of these 
peculiarities, the fermionic propagator has the expected l/$f pole on the 
domain wall. But there may be a problem with the phase of the fermionic 
determinant at the non-perturbative level. 
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Numerous attempts [1] in the past to define chiral gauge theories using lattice 
regularization have been unsuccessful because of the doubling problem [2]. Recently, 
Kaplan [3] proposed that a lattice chiral gauge theory in 2n dimensions could arise 
as the low energy limit of a massive Dirac theory in In + 1 dimensions, provided the 
fermion mass takes the form of a domain wall 

m(s) = m sign(s) , (1) 

where s denotes the extra coordinate. Besides massive (2n+l)- dimensional excitations, 
the hamiltonian has a spectrum of bound states [3-5] 

i/>(x,s;p) = e - mlsl+tp - x u{p) . (2) 

Here denote the 2n-dimensional coordinates and u{p) is a positive chirality spinor 
satisfying $u(p) = 0. Eq. (2) describes a massless chiral fermion that lives on the 
domain wall. 

A naive discretization of the continuum hamiltonian does not give rise to chiral 
fermions because the lattice zero mode has the expected right handed doubler. How- 
ever, Kaplan [3] and Jansen and Schmaltz [6] showed that by adding a Wilson term 
one can obtain a chiral spectrum for a certain range of the ratio m/r where r is the 
Wilson parameter. For < \m/r\ < 2 the free lattice hamiltonian describes one Weyl 
fermion near the origin of the Brillouin zone. For higher \m/r\ the chiral fermions 
move to the vicinity of other corners of the Brillouin zone. 

Thanks to current conservation which is always true in odd dimensions, the 2n- 
dimensional anomaly can be written as d^j^ = — $2n+]j2n+i- Thus, the anomaly arises 
as a non-zero flux away from the domain wall. The flux has been calculated [7] for 
general m/r in the presence of an external, slowly varying gauge field. It is a non- 
trivial result that the flux takes only integral values, and it changes discontinuously as 
a function of m/r in agreement with the anomaly anticipated from the lattice chiral 
spectrum. 

In this paper we consider the Dirac operator of the euclidean theory. Under- 
standing the properties of the Dirac operator is important for obvious reasons. The 
euclidean partition function is better defined mathematically than the Minkowskian 
path integral, and numerical simulations rely heavily on the euclidean metric for 
convergence. 

Our discussion is limited to the Dirac operator of the free theory, its only non- 
trivial feature being the domain wall shape of the fermion mass. We first discuss the 
properties of the continuum Dirac operator. Surprisingly, we find that its spectrum is 
drastically deferent from the spectrum of the hamiltonian. The most important result 
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is that there exist no bound states with non-zero momentum! Another peculiarity is 
that the zero momentum spectrum contains an infinite number of bound states. For 
general momentum, the spectrum consists of eigenstates which are superposition of 
plane waves in one half of the (2n+l)-dimensional space and which decay exponentially 
in the other half. Only these eigenstates contribute to physical observables because 
Pn = is an isolated discontinuity point of the spectrum. 

The peculiar behaviour described above is possible because the euclidean Dirac 
operator is not hermitian. For hermitian operators such as the hamiltonian, the exis- 
tence of a zero momentum bound state separated by an energy gap from the massive 
modes implies that finite momentum bound states must also exist. However, the 
variational argument needed to prove the above statement fails when the eigenvalues 
are complex. 

On the lattice we discuss in detail the case r — 1. Taking a finite lattice with 
periodic boundary condition in the s direction we find yet another surprise. For 
general m, even the p M = spectrum contains no (approximate) zero mode. Most of 
the eigenstates do look like bound states, but the eigenvalues are typically all 0(1). 
The only exception is when m takes the special value m = In this case one has 
a bound state spectrum with no energy gap, but the small p dispersion relation is 
quadratic rather than linear. (The exact form of the dispersion relation is non-local). 
The case r ^ 1 is technically more complicated and we do not discuss it in detail. We 
expect that the qualitative properties of the spectrum should be the same as in the 
r = 1 case. 

The eigenstates of the euclidean Dirac operator do not have simple physical in- 
terpretation. What enters the calculation of physical correlation functions is the 
fermionic propagator and the fermionic determinant. The continuum fermionic prop- 
agator can be constructed by considering the operator D F D F . Denoting by Q the 
propagator of DpD F , the fermionic propagator Qp is given by Qp = D F Q. 

For p^ = 0, DpDp has a negative chirality zero mode. Since DpD F is hermitian, 
the zero momentum bound state implies the existence of a bound state spectrum for 
p^ 7^ as well, which, in turn, gives rise to a 1/p 2 pole in Q on the domain wall. 
Consequently, Qp has the desired chiral 1/p 1 pole on the domain wall in spite of the 
absence of finite momentum bound states in the spectrum of Dpi 

On a finite lattice, Dp is a matrix and DpD F is defined through matrix multi- 
plication. Because of the Wilson term, the lattice D F D F may contain up to fourth 
lattice derivatives, and so arguments based on the operator D F D F are more delicate. 
In the case r = 1 this difficulty does not arise. The lattice DpD F contains only second 
derivatives, and its only unusual feature is that the normalization of the kinetic term 
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is different in the regions s > and s < 0. Thus, the argument goes through and, at 
least in the r = 1 case, we can prove that the lattice propagator has a chiral 1/p 1 pole 
on the domain wall when m is in the allowed range [6]. 

A similar argument applies to the modulus of the fermionic determinant. It, too, 
can be constructed from DpD F because |Det{Z?ir}| = Det^ {DpD F }. When gauge 
fields are introduced, one can expect that the response of fermionic propagator and 
the modulus of the fermionic determinant will be dominated by the chiral mode of 



The object which must be obtained directly from Dp is the phase of the fermionic 
determinant. The peculiar spectrum of Dp may indicate an unexpected behaviour of 
the phase of Det{D F }. Naively, one can imagine several scenarios according to which 
the unusual properties of the euclidean spectrum may prevent the existence of a chiral 
continuum limit. Finding what really happens requires a detailed understanding of 
the response of the phase of Det{D F } to gauge fields. We comment that the Chern- 
Simons action found in ref. [7], which reproduces the usual anomaly cancellation 
condition, arise in weak coupling perturbation theory. However, we expect that the 
continuum limit should depend on the dynamical properties of the five dimensional 
gauge theory, and hence that non-perturbative contributions should play an important 
role. 

We begin our analysis with the continuum Dirac operator 



For simplicity we work in (4 + 1) dimensions. All gamma matrices are hermitian and 
repeated indices are summed from one to four. Most of the discussion generalizes 
trivially to other dimensions. The dependence on dimensionality will be explicitly 
given whenever necessary. 

Assuming a plane wave solution in the first four coordinates we obtain an effective 
one dimensional problem 



DpDp. 



D F = il^d^ + ry 5 <9 s + im(s) . 



(3) 



Dp = ry 5 <9 5 + im(s) + f 1 . 



(4) 



cece We consider the eigenvalue equation 



Dpij) = Xijj . 



(5) 



The general solution can be written as 



tp = U(s;X,p)u(X,p,a) , 



(6) 
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where U(s; X,p) is a four by four matrix and u(X,p, a) is an s-independent spinor. a 
is a spin index which we will usually suppress. Integrating eq. (5) we find 



cosh(s/t + ) + k + 1 75(?(^ — A) — m)sinh(s/t + ) , s > , 
cosh(s/t_) + KZ l ^i^{i{i> — A) + m)sinh(s/t_) , s < , 



where 

k 2 ± = (iX ± m) 2 + p 2 . (8) 

For definiteness, we pick the branch were k± — > iX ± m for p 2 — > 0. Notice that U(s) 
is continuous at s = and has a discontinuous derivative there as required by eq. (4). 
One can obtain an analytic dependence on s by smoothing the discontinuity of the 
mass function at s = 0. However, the qualitative properties of the euclidean spectrum 
depend only on the asymptotic behaviour of m(s) at ±oo, and they are completely 
insensitive to the details of the transition region. 

For general values of A, the solution U(s) contains exponentially increasing and 
exponentially decreasing pieces. We begin our investigation by looking for special 
values of A for which the solution is a superposition of plane waves in some region of 
the s-axis. Demanding for example Ken + = we find that A must take the special 
form A = £ + im where the real parameter £ satisfies £ 2 > p 2 . In the region s > 0, 
U(s) is then a superposition of plane waves with k = ±(£ 2 — p 2 )^ . To obtain U(s) for 
s < we need k_. Assuming £ 2 , p 2 m 2 we find ~ —2m — p 2 /4m + i£. Hence, 
the solution has no plane wave component in the region s < 0. For given p^ and A 
there are four linearly independent solutions, of which two increase exponentially and 
two decrease exponentially for s < 0. These solutions are obtained by an appropriate 
choice of the constant spinor u(X,p), see below. 

The spectrum is obtained by first taking s to lie in the finite interval — L < s < L 
with some local boundary conditions at s = ±L. We obtain the spectrum below 
assuming periodic boundary conditions. (Other choices of the boundary conditions 
lead to similar conclusions). In the limit L — > oo, only solutions which decrease 
exponentially for s < remain in the spectrum. Of course, there is another family of 
solutions for which the role of s > and s < is interchanged. 

Let us now examine the solution of eq. (5) for |ImA| < m. In this case one has 
Ren + > and Re/t_ < 0. If no special choice of u(X,p) is made, the solution will 
contain exponentially increasing and exponentially decreasing pieces for both s > 
and s < 0. To obtain only exponentially decreasing (respectively increasing) solution 
for s > we demand 

A>u(X,p)=0, (9) 
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where 



> 1 i(^-X)-m 



Notice that are orthogonal projection operators. Similarly, taking into account 
that ReK_ < 0, an exponentially decreasing (respectively increasing) solution for 
s < is obtained by demanding 

A<u(A,p)=0, (11) 

where 

A < = i ±T5 3t^n (12) 

A normalizable eigenstate is obtained by demading that the solution decrease 
exponentially on both sides of the domain wall, i.e. u(X,p) must satisfy A>u(\,p) = 
A<u(\,p) = simultaneously. Using eqs. (10) and (12) we find that a non-trivial 
solution exists if and only if the matrix (k+ — K- — 2m7 5 ) has a zero eigenvalue. For 
this to happen, we must have = 2m, which, in turn, is true only when p^ = 0. 

In this case (k + — k_ — 2m7 5 ) becomes a projection operator and the bound state 
has a definite chirality. 

We therefore arrive at the conclusion that euclidean bound states exist only for 
Pfi = 0. When p^ ^ 0, one can obtain solutions that decrease exponentially on one 
side of the domain wall but not on both. The technical reason for this behaviour 
is that, as can be seen from eq. (4), in euclidean space f couples the positive and 
negative chiralities. (By contrast, the kinetic term of the hamiltonian a -p commutes 
with 75). The homogeneous solution with negative chirality increase exponentially on 
both sides of the domain wall, and for finite p^ one can eliminate the exponentially 
growing component only on one side. 

For p^ = the positive and negative chiralities decouple, and U (s) simplifies to 

U(s) = ^(1 + 75 )e- m l s HA S + 1(1 - l5 )e m ^ +iXs (13) 

We readily see that for p^ = there is a positive chirality bound state for all A 
provided |ImA| < m. 

Our next step is to obtain the spectrum by imposing periodic boundary conditions 
in the finite interval — L < s < L, later taking the limit L — > 00. Imposing periodic 
boundary condition amounts to introducing an anti-domain wall at s — L. By the 
same argument as before, bound states which live on the anti-domain wall exist only 
for p^ = 0. These bound states correspond to the negative chirality part of eq. (13). 

Eigenstates which satisfy periodic boundary conditions occur when the matrix 
equation 

(U(s = L) - U(s = -L))v = 0, (14) 
6 



has a non-trivial solution. The zero eigenvectors are than identified with the spinors 
u(X,p, a). We must therefore look for the zeros of det (U(L) — U(—L)). After some 
algebra we find in 2n + 1 dimensions 

det (U(L) - U(-L)) = A 2 ™ -1 , (15) 

where 

1 p 2 — w? — a 2 

-A = 1 — cosh(L/t + )cosh(— Lk_) H sinh (Lk + ) sinh (— LkJ) . (16) 

2 tx— 

Consider first the special case p^ = 0. In this case further simplification occurs 
and we obtain 

^A = 1 - cosh(2iAL) . (17) 

The zeros occur for A& = rck/L. Thus, imposing the periodic boundary conditions 
eliminates all bound states with ImA ^ from the spectrum. However, we are still 
left with an infinite number of bound states, which becomes a continuum of bound 
states where A is any real number in the limit L — > oo. 

The existence of infinitely many bound states for p^ = provides another expla- 
nation for the absence of finite momentum bound states. Had every p^ = bound 
state been the end point of a continuous four dimensional spectrum, we would have 
found that the number of fields that live on the (anti)-domain wall is infinite! We 
comment that, the fact that the eigenvalues lie on the real axis is a consequence of the 
symmetry between the two sides of the domain wall. Had \m(s > 0)| been different 
from \m(s < 0)|, we would have found that the bound states' eigenvalues move away 
from the real axis. We will see that a similar phenomenon takes place on the lattice. 

We next consider A for p^ ^ 0. We will show that the zeros of A occur for |ImA| ~ 
m + 0(1/L). Let us first assume that |ImA| — m is large in units of 1/L. In this case 
A is dominated by the exponentially increasing part of the cosh and sinh functions. 
Therefore, a necessary condition for having A = is that the difference between 
the coefficient of the cosh term and the coefficient of the sinh term be exponentially 
small. Explicitly, we must have k\ k 2 _ = (p 2 — m 2 — A 2 ) 2 + 0(e~ mL ). This condition 
is satisfied only if p/m m e~ mL . 

In the hamiltonian picture, boundary effects on the chiral mode are 0(e~ mL ). 
Thus, for exponentially small momentum we expect a significant distortion of the 
chiral mode and so p/m e~ mL is of no physical interest. 

Assuming p/m^> e~ mL , we conclude that A has no zeros if |ImA| — m is large. The 
zeros of A must therefore be found in the vicinity of |ImA| = m. For definiteness let 
us assume ImA ~ m. A convenient parametrization is A = ^+i{m — kf3/£ i L). (We take 
k = +(i 2 —p 2 )^ 2 and £ 2 > p 2 as before). This implies k + ik + fi/L. Now, cosh(L/t + ) 
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and sinh(LK + ) are 0(1), whereas cosh(— LkJ) ~ sinh(— Lk) ps exp(— Lk_). (Notice 
that Re/t_ is negative). Hence 

/ 2 _ m 2 _ ^2 \ 

A ~ exp(— Lk_) cosh(L/t + ) H sinh(L/t + ) . (18) 

V J 

Up to exponentially small corrections, the zeros of A occur for 

tgh(L*; + )= 2 K+K ~ (19) 
p 2 — m 2 — \ 2 

We can find the zeros explicitly in the limit p 2 ,C, 2 <C m 2 . We first notice that 

1/r . sinh(28) + isin( 2 Lk) , nn . 

whereas 

2 ^' x2 = ^ + 0(l/m). (21) 
p 2 — m 2 — X 2 £ 

Since fc/^ is real, the zeros occur when sin(2/cL) = 0. We than determine (3 from the 
equation 

sinh(2/3) _A; 



cosh(2/3) ± 1 f 

The finite L eigenstates have the following shape. The spinor u(X,p, a) is made 
almost entirely out of the projection, plus an exponentially small contribution 
from the Af projection. When we go from s = into the negative s region, the wave 
function first decreases exponentially. Around s = —L/2 the other solution takes 
over, and the wave function starts growing exponentially, until at s = —L it becomes 
0(1) again. In the limit L — > oo, (ImA — m) vanishes, the exponentially growing 
region recedes to — oo and we obtain the infinite space spectrum described earlier. 

We now turn to the lattice Dirac operator. We take a lattice of 2N sites in the s 
direction with periodic boundary conditions, i.e. the points s = N and s = —N are 
identified. We set m(s) — m for s — 1, . . . , AT and m(s) = —m for s = —N + 1, . . . , 0. 
(Alternatively, we could take m(s) = at the sites s = and s = N. For large A^ 
this modification has a negligibly small effect). 

As in the continuum case we assume a plane wave solution in the first 2n coordi- 
nates. The effective Dirac operator in the s direction becomes an 8A^ x 8A^ matrix. 
With lattice spacing a — 1, the eigenvalue equation reads 

^75 (U(s + 1) - U(s - 1)) + r - (U(s + 1) + U(s - 1) - 2U(s)) 

+ h^l^m{P^-rF{p)+m{s)+tXJU(s) = 0, (23) 



2n 

F(p) = £(l-cos(^)). (24) 
(t=i 

Here U (s) is a four component spinor. As before A is the euclidean eigenvalue. 

We consider first the case r — 1. Denoting by U±(s) the positive (negative) 
chirality components, Eq. (23) simplifies to 

2n 

U + (s + l) = (l + F{p)-i\-m{s))U + {s)-iJ2v»sm(P»)U-{s), (25a) 

2n 

U-(s-l) = (l + F( P )-i\-m(s))U-(s)-iJ2°^HPn)U + (s). (25b) 

We introduce transfer matrices defined by 

U{s + 1) = T + U(s), s>0, (26a) 
U(s-l) = T~U(s), s<0. (26b) 

To obtain e.g. T + we use eq. (25b) to express LL(s+l) in terms of LL(s) and U + (s+l), 
and eq. (25a) to express U + (s + l) in terms of U + (s) and £/_(s). We find 



7 1 + F(p) -iX-m sin(p M ) X 

v S^^sin^) ^gg-^ 



(27) 



where 



2n 



G(p) = ^sin 2 (p M ). (28) 

Notice that det(T + ) = 1. To obtain T~, invert T + and replace m by — m. 

We will analyze the spectrum in the vicinity of the origin of the Brillouin zone. 
We assume < m < 2 so that on an infinite lattice the zero mode [3] indeed occurs 
at = 0. For larger values of m the zero modes [6] occur at other corners of the 
Brillouin zone, but the analysis remains essentially the same. 

We begin with the = case. The positive and negative chiralities decouple 
and we consider the positive chirality eigenvectors. We will show that the positive 
chirality zero mode [3] of the infinite lattice disappears on a finite lattice (no matter 
how large) unless we fine tune m to a special value. In the positive chirality channel, 
the periodicity requirement reads 

(l-i\-m) N = (l-i\ + m)- N . (29) 

The eigenvalues which solve eq. (29) are 

iX± = l±(m 2 + e^f , k = 0,...,N-l. (30) 
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The eigenvectors are given by U + (s + l) = (l — i\ — m)U + (s) for s > and U + (s — l) = 
(1 — iX + m)~ 1 U + (s) for s < 0. Notice that most of the eigenvectors look like bound 
states which live on one of the domain walls. 

A zero mode occurs only if N is even and m = y/2. For arbitrary to, the eigen- 
values are all 0(1). This has the following explanation. On an infinite lattice, a 
euclidean bound state exists for any A which belongs to the intersection of the open 
sets defined by |1 — iX — m\ < 1 and |1 — iX + m\ > 1. (For to in the allowed range, 
the intersection always contains a neighbourhood of the origin in the complex plane). 
This region is the analog of the stripe |ImA| < m in the continuum case. The finite 
lattice's boundary conditions pick a finite number of points out of this region. The 
Wilson term breaks the symmetry between the two sides of the domain wall, and so, 
unlike the continuum case, the eigenvalues do not have any special reality property. 

Since the lattice Dirac operator is a finite matrix, its eigenvalues and eigenvectors 
solve algebraic equations and so they must depend continuously on p^. The magnitude 
of the deviation from the p M = behaviour is controlled by the parameter \p^N\. For 
<C N^ 1 the distortion will be small, whereas for p^ > N^ 1 we expect that the 
shape of the eigenvectors will be qualitatively different from the p^ = case. 

We now want to find whether the lattice supports a spectrum of bound states for 
N^ 1 < p M < 1. The analysis parallels the continuum case (although the results are 
different) and we only outline it here. The transfer matrix T + has two eigenvalues 
which satisfy rji = l/rj}- (In 2n + 1 dimensions the multipliticity of each eigenvalue 
is 2 n ~ 1 ). Without loss of generality we assume that l^+l < 1 . A similar statement 
applies to T~. 

We define projection operators into the subspaces that belong to each eigenvalue 
separately for T + and T - . In order to have a bound state, the rji and rji subspaces 
must have a non-trivial intersection. Calculating to leading order in p^ and A we find 
the necessary condition 

2-m 2 + 2p 2 -2iX^0. (31) 

Eq. (31) holds only if m = y/2 and iX xs p 2 . In order to verify that the bound states 
belong to the spectrum we must check that they satisfy periodic boundary condition. 
It is easy to verify that, provided 

i\ = l + F(p)-(l-G(p))* , (32) 

one has T + = — T~. Thus, taking iV to be even we conclude that the bound states 
belong to the spectrum and that the exact dispersion relation is given by eq. (32). 
We see that, in the special case where the lattice spectrum has an exact zero mode, 
this zero mode is the end point of a bound state spectrum. 
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The bound state spectrum found above describes a lattice excitation that lives 
on the domain wall. But it is clear that this mode does not describe a relativistic 
fermion. What the euclidean mode (32) has in common with the chiral mode of 
the hamiltonian is (a) the absence of an energy gap (b) the existence of two states 
per given momentum and (c) the absence of other zeros except the one at p^ = 0. 
But there are also important differences. The euclidean mode has a definite chirality 
only at p^ = 0. The absence of an energy gap depends on fine tuning the mass 
parameter, whereas the masslessness of the chiral mode of the hamiltonian arise 
because of topological reasons. Finally, the dispersion relation (32) is non-local, and 
for small p^ it is quadratic rather than linear. 

The mass function (1) is of particularly simple shape, and one could consider the 
effect of a more general ansatz. Since the fifth direction is not physical, in the most 
general case the mass function may depend on chirality and may contain non-local 
couplings in the s direction. The only absolute requirement is that m(s) respect 
locality and spacetime symmetries in four dimensions. Remarkably, the behaviour 
exhibited by the mode (32) seems to persists in the most general case in euclidean 
space. 

In trying to understand this behaviour, we notice that properties (a-c) above 
can be satisfied without even introducing a fifth dimension. All That is needed is to 
take the four dimensional massless Dirac operator and to add a "chiral" mass term 
|(1 — 75)772. The reader can easily check that the small p^ dispersion relation is then 
quadratic instead of linear. Obviously, in this case the no-go theorems apply, and 
the fact that the small p^ dispersion relation is quadratic is crucial in order to avoid 
additional zeros inside the Brillouin zone. 

In the present context we have a five dimensional problem. However, since we 
are dealing with free fermions, picking a particular mode automatically defines an 
effective (in general non-local) four dimensional theory. Concentrating on the small 
p^ behaviour, whenever we demand that this mode satisfy properties (a) and (b) 
above, we end up with an effective Dirac operator that has the "chiral" mass term, 
with the resulting quadratic dispersion relation. 

We comment that the lattice hamiltonian evades the no-go theorems because of 
the possibility of level crossing. As in four dimensions, the spectrum must be periodic 
across the Brillouin zone, but when we keep track of the flow of a particular eigenvalue 
as a function of p M , we find ourselves at p = 2n on a different branch from the one we 
started with. We have already seen that euclidean bound states have little to do with 
the bound states of the hamiltonian. It will be interesting to check whether in the 
present context one can prove a no-go theorem only for the euclidean Dirac operator. 
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We finally comment that the case r ^ 1, while being technically more complicated, 
readily contains all the ingredients that make the euclidean spectrum so different 
from one's naive expectation. In particular, we expect that most of the small p M 
eigenvectors will look like bound states, but that a massless mode will occur only 
incidentally. Furthermore, whenever such a mode should occur, it is bound to suffer 
from the same problem as discussed above, i.e. to have a quadratic instead of a linear 
dispersion relation. 

After this lengthy discussion of the peculiarities of the euclidean spectrum, the 
demonstration that the fermionic propagator has the desired l/$f pole on the (anti)- 
domain wall is remarkably simple. We consider first the continuum case. One has 

D F D F = -d 2 - + m 2 + 2m l5 5{s) . (33) 

In an infinite space D F D F has a bound state spectrum 

ij(x,s;p) = e- mlsl+tp - x u(p) , (34) 

where u(p) is an arbitrary negative chirality spinor, and 

D F D F i[>(x,s;p) = p 2 ip(x,s;p) . (35) 

The dispersion relation is therefore A 2 = p 2 . It is strictly local, reflecting the separa- 
tion of variables which characterizes the mode (34). The contribution of the mode (34) 
to the propagator of D F D F is 

i „— m\s\ — m\s'\ 

g chiral (s, s';p) = -(l- 7s )- . (36) 

z p 

Since tt (1 - 75 )e- m l s l is annihilated by D F , the contribution to the fermionic propa- 
gator is 

gf ral (s,s';p) = D F g chira \s,s'-p) 

i -m\s\-m\s'\ 

= 2 {l+ ^^^- < 37 > 

Clearly, eq. (37) describes the propagation of a chiral fermion on the domain wall. 

When we impose periodic boundary conditions, the propagator develops another 
pole which describes a Weyl fermion with the opposite chirality that lives on the 
anti-domain wall. The crucial point is that, using the hermiticity of D F D F we can 
rigorously show that the chiral poles are modified at most by 0{er mL ) corrections. 
In fact, in the present case there is no modification at all because the mode (34) is 
symmetric under s — > — s. 
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As we have explained in the introduction, similar conclusions apply to the lattice 
propagator in the case r = 1. On the lattice the boundary effects are 0(e~ N ). They 
do not vanish identically because of the absence of an s — > —s symmetry. We have 
not analyzed in detail the case r ^ 1. In any event, since the lattice effects tend to 
zero in the limit of small m and p^, we expect that the continuum behaviour should 
persist for some finite range of m, if not for the entire interval found in ref. [6]. 

We conclude with a summary of our main results. Any object which can be 
obtained from the hermitian operator D F D F is "well behaved". In particular, the 
fermionic propagator has the expected chiral poles at least for some range of m and 
r values. Since perturbation theory is defined as the sum of one loop diagrams, 
we expect that weak coupling perturbation theory on the lattice should lead to the 
same conclusions as found when perturbation theory is defined through dimensional 
regularization. Thus, our conclusions support the results of ref. [7]. 

At the non-perturbative level we still expect the modulus of the fermionic deter- 
minant to be well behaved. But the phase of the fermionic determinant may exhibit 
a peculiar behaviour because it must be constructed directly from the spectrum of 
D F itself. 

To make our point more clear, consider the construction of the fermionic propaga- 
tor directly from the spectrum of Dp. We remind the reader that, being complex, Dp 
has "right eigenstates" and "left eigenstates" . So far we have discussed only the right 
eigenstates, which are defined by the eigenvalue equation (5). The left eigenstates are 
obtained when we take Dp to act to the left, i.e. 

XD F = \x- (38) 

(On the lattice, the left action of Dp is represented by the transposed matrix). The 
left and right eigenstates are mutually orthogonal 

[ a /3 ! \ — Aj . 

The indices a and f3 count the multipliticity of each eigenvalue. The propagator has 
the generic form 

^_^ A a (x,s)xgx',s') ^ (4o) 

i,a A * 

As we have explained earlier, the pole in the propagator should persist own to 
p e~ mL . This means that in the limit of small p M , the norm of Qp can be very large. 
In the continuum case, the propagator is constructed from the balk eigenstates of the 
two half-spaces. The interesting region is e~ mL <C p/m <C I. Since |ImA| rs to, the 
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denominator on the r.h.s. of eq. (40) is never small. This observation means that, in 
order to build the pole singularity, balk eigenstate up to ReA ~ m 2 /p must contribute 
coherently to the fermionic propagator in the vicinity of the domain walls. Thus, the 
ultraviolet and the infrared behaviour of the euclidean spectrum are closely related. 

On a finite lattice this mechanism cannot work, because the number of eigenstates 
per given momentum is just 8iV, whereas 1/p can be as large as e N . How does 
the r.h.s. of eq. (40) manage to be exponentially large? The crucial observation 
is that the mutual orthogonality condition (39) does not imply that ip and x are 
normalized. Since the terms that contribute to the chiral pole are of the form xip^i 
their magnitude is not constrained by eq. (39). In fact, on the lattice the off-diagonal 
terms Xl^ are exponentially large whenever the eigenvectors look like bound states. 
As we have shown, for p^ <C iV -1 most of the eigenvectors indeed look like bound 
states. Thus, the pole on each domain wall receives contributions from approximately 
half of the eigenstates. 

The physical picture that emerges from our discussion is the following. In eu- 
clidean space, the chiral propagation on the domain wall is a coherent effect of the 
entire five dimensional spectrum. No single mode of the lattice Dirac operator is 
capable of producing the chiral pole. 

Turning to the interacting theory we comment that, roughly speaking, what is 
needed in order to have a chiral continuum limit is that the fermionic determinant 
be factorizable into four pieces. Two of them should depend on the four dimensional 
gauge fields in the vicinity of the (anti)-domain wall, while the other two should 
depend on the gauge degrees of freedom in the two balk half-spaces. If this scenario 
works, we can expect that a continuum limit should exists and that it should describes 
two chiral theories (one on the domain wall and one on the anti-domain wall) which 
no nothing of each other. What might go wrong as a consequence of the peculiarities 
of the euclidean spectrum is, for example, that non-local correlations will develop in 
the s direction. Another interesting possibility is that a chiral continuum limit does 
exist, but that its dynamics is qualitatively different from QCD. Finding what really 
happens must await further investigations. 
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